Charge pumping in a quantum wire driven by a series of local time-periodic potentials 
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We develop a method to calculate electronic transport properties through a mesoscopic scattering 
region in the presence of a series of time-periodic potentials. Using the method, the quantum charge 
pumping driven by time-periodic potentials is studied. Jumps in the pumped current are observed 
at the peak positions of the Wigner delay time. Our main results in both the weak pumping 
and strong pumping regimes are consistent with experimental results. More interestingly, we also 
observed the nonzero pumping at the phase difference <f> — and addressed its relevance to the 
experimental result. 
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A parametric electron HUaap has attracted considerable 
attention in recent yearsJj~£3 It is a device that generates 
a dc current at zecp bias by cyclic deformations of sys- 
tem parameters.tTcl The quantum pumping mechanism 
was originally proposed by Thouless,D who studied the 
integrated particle current produced by a slow periodic 
variation of the potential, and showed that in a finite 
torus the integral of the current over a period can vary 
continuously, but it must have an integer value in an 
infinite periodic system with full bands. Such quantized 
charge transport was proposed to become an electric cur- 
rent standard.u 

Quite recently, the charge pumping was observed ex- 
perimentally.!] For technical reasons, instead of measur- 
ing charge currents, the pumped dc voltage Vdot is mea- 
sured in a quantum dot where two gates with oscillat- 
ing voltages control the deformation of the shape of the 
dot. For weak pumping, the observed charge pumping 
has a sinusoidal dependence on the phase difference <f> 
between the two shape-distorting ac voltages applied to 
the gates, and is proportional to the square of pumping 
strength V. For strong pumping, the pumped current 
deviates from the square dependence on V and becomes 
nonsinusoidal, being always antisymmetric about <f> = n. 
The charge pumping may have a close relation to the 
adiabatic Berry's phase since the evolution of the sys- 
tem is cyclic and is controlled by several system param- 
eters, referred to as the parametric pumping. Based on 
this understanding, the total charge pumped per cycle 
is proportional to the area enclosed by the path in the 
parameter space, and nonjzero pumping current requires 
at least two parameters. Era The pumped charge drived 
by two parameters should be zero if two parameters are 
in phase (0 = 0)) since the area enclosed by the path is 
zero. However, it is in contradiction with the observed 



current I((f> = 0) ^= O.i One of possible mechanisms of 
nonzero currents for cf> — is photovoltaic effects intro- 
duced in Ref.cl, where a surprising result, nonzero dc cur- 
rent generated by a single pumping gate voltage, is also 
reported. The general physics of a quantum PWBP has 
been the subject of several theoretical analyses. Em Zhou 
et al demonstrated that at low temperatures both the 
magnitude and the sign of the pumped charges are sam- 
ple specific quantities, and the typical value in disordered 
(chaotic) systems turns out to be determined by quan- 
tum interference effects. Another general expression for 
the average transmitted charge current was derived by 
BrouwerH under the adiabatic condition and based on 
the time-dependent S-matrix method,Ll which appears to 
be quite successful for (adiabatic) weak pumping. Adia- 
baticity here means that the oscillating period r of 
system is much larger than the Wigner delay time t w 
Note that the adiabatic condition does not simply im- 
ply that the pumping strength V should be very small. 
In fact, the adiabatic condition requires that r must be 
larger as V increases. On the other hand, the pumping 
was not weak in the experiments!] The main purpose of 
the paper is to develop a theory, which is also applicable 
in the case of strong pumping. By using the Floquet the- 
orem, the. photon-assisted transport has been taken into 
accounttil. We calculate the pumped current through a 
mesoscopic region in the presence of time-periodic poten- 
tials. Our main results in the weak pumping regime, as 
well as those in the strong pumping regime are consistent 
with the experiment reported in Rcf.El. 

Consider electrons transmitting through a one- 
dimensional scattering region ranging from xo to xq + S. 
The potential is given by 



V(x,t) 



0, xo < 0, x > xo + S 

V s (x,t), Xq < X < Xq + S 



with V s (x, t) = Vo + V s cos(u>t 
equation can be written as 



(1) 

). The Schrodinger 
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n 2 d 2 ^(x,t) 



2m* 



dx 2 



+ V(x,t)9(x,t), (2) 
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with m* as the electron effective mass. Eauation g can 
be solved by using the Floquet theorem .113 By setting 
ty Fl (x,i) = e- iEFlt / h ip(x,t), where E F i is the Floquet 
eigenenergy and ip(x,t) is a periodic function ip(x,t) — 
ip(x,t + r) with period r = 2ir/u>, the Schrodinger equa- 
tion takes the form 

. . . ft 2 d 2 tb(x,t) ^dip(x,t) T ,, , ,, 
E Fl ^(x,t) = - — 2* > - t h-!?yJ-+v(x,t)il>(x,t). 



dt 



Substituting ip(x,t) — g(x)f(t), we have two separated 
equations with an introduced constant E, 



h 2 d 2 

2m* dx 2 



g(x) + V g(x) = Eg(x), 



(3) 



iti 



df(t) 

at 



V s cos(cot + fa)f(t) = (E- E Fl )f(t). (4) 



Integrating Eq.(||) gives 



iV a si n <t> s i(E-E Fl )t 

f(t) = e «<- 



— ; 

raw ' 



\ (5) 



where J„(a;) is the Bessel function of the first kind of 
order v. Since f(t) is periodic in time with period r, it 
follows from Eq.(^|) that E — Epi = muj with m as an 
integer. The equation for g(x) has a solution 

g(x) = e ±ik ™ x , (k s m ) 2 = 2m*(E Fl +mhuj-V )/h 2 . (6) 
Thus ip{x,t) becomes 



ip m (x, t) = e *(Vi/»w)«in*.±a&» J2 F n 



\ (7) 



with F„_ m = ea;p[-i(n - m)fa] J n - m (V s /huj). 

We consider an incoming wave from the left with the 
energy Eq = % 2 k^/2m* , then the outgoing waves should 
be divided into different modes E n , which satisfies E n = 
E +n?iuj with n = 0, ±1, ±2, ■ • -. The propagating modes 
mean that E n > 0, while the evanescent modes mean that 
E n < 0. The latter exists only in the neighborhood of the 
oscillating barrier and do not propagate. Denote k n — 
\/2m* E n /h, the solutions of the Schrodinger equation 
can be written as 

oo 

= E {^ lknX + K^ lknX )^ lEnt,n , {x < X ) 

71— — OO 

oc 

*. = e-^ £ (a m e ik - x + b m e-^ x )F n ^ m e- m ^, 

m,n— — oo 

( Xq < X < Xq + 5) 

oo 

* r = E (B % n e~ iknX + B°e iknX )e- iEnt / h , (x > x + 5) 

where A l n and B % n are the probability amplitudes of the 
incoming waves from the left and right, respectively, 
while A° n and B° are those of the outgoing waves. We can 



characterize the barrier by a scattering matrix S which is 
a matrix connecting the incoming and outgoing channels 



A° 
B° 



S 



B' 



where the S matrix can be derived by the matching con- 
ditions for the wave function ^>(x,t) and its derivative 
d x ^{x,t) at x = lEor-apd x — xq + S. After eliminating 
a m and b m , we havefij 



S = 



(8) 



LpTLp 1 , R^ 



LrRLr, T- 



where T«_ - ^ ^ R , - ^ R ^ R ,^^-^ R 
and = L^RLl. Here the left (right) arrow indicates 
incoming waves from right (left), the matrices Ll and Lr 
are defined as {L L ) mn = exp[ik n x ]S mn and (L R ) mn = 
exp[ik n (xQ + S)]5 mn . T and R are given by 



where 



T = (C^ 1 D 1 + C 2 ~ 1 D 2 )/2, 
R = (C 1 7 1 D 1 - C 2 ~ 1 D 2 )/2, 



C x = (L s - I)K S F^ - (L s + T)F^K, 
D x =-{L s - I)K S F^ - {L s + I)F^K, 
C 2 = (L 3 + I)K S F^ - (L s - I)F*K, 
D 2 = (L s + I)K S F^ + (L s - TjF^K, 



(9) 
(10) 



with the matrices (L s ) mn ^ = exp[ik s n S\5 mn , (K s ) mn = 
k^Smn, K mn = k n S mn , I as the unit matrix and F> 
as the Hermitian conjugate of F. The electronic trans- 
port properties of the scattering region may be obtained 
straightforward from Eq.(^j). 

The above method may be generalized to I time- 
periodic barriers described by 



V(x,t) 



0, X < 0, X > CLl 

Vi(x,t), < x < ai 
V 2 (x, t), a\ < x < a 2 

. Vi(x,t), a ; _i < x < ai 



(11) 



where Vi(x,t) = V? + Vicos(u;it + fa), V 2 (x,t) = V 2 ° + 
V 2 cos{uo 2 t + fa), • • • and Vi(x, t) = V t Q + Vicos(u>it + fa). 
This potential may be more a realistic model for experi- 
ments. Obviously the transport properties for each bar- 
rier can be characterized by an S matrix given by 



S a 



R'l 



rpo 

R" 



where a = 1, 2, • • • , I, T°, TZ, R°U and R?_ can be de- 
rived by the same method presented above. Now the 
propagating mode E n should be replaced by 



2 



E( nj )=E + n ^r (i = (12) 

rij — — oo 

The associated transfer-matrix M a for the ath barrier 
may be derived directly from the S a matrix 

The total transfer-matrix M* for all those barriers is de- 
termined by 

where = 1,2) are the partitioned matrices with 

the same size as Tf;. The total scattering matrix S* can 
be derived from M l as 

In each cycle a net charge current may pass through 
the scattering region in the direction determined from the 
detailed form of 5*' matrix. We define a net transmission 
coefficient (for an incoming wave in mode Eq — E) by 

T «*= E ^^(l^o(^o)| 2 -|Tl,„ 30 (i?o)| 2 ). 

E( nj )>0 

The average net current per period r (for Eq) through 
those barriers is j(Eo) = T net (Eo). If the system is con- 
nected through two ideal leads to two electron reservoirs 
with the same chemical potentiaLu^the average pumped 
current per period r is given bytLj'tj 

/>oc 

I(p) = e / dEg(E)f(E - fx)T net (E), (13) 
Jo 

where g(E) — ^J2m* /E/h is the density of electrons con- 
tributing to the current in one direction, and f(E — /i) 
is the Fermi-Dirac distribution. At zero temperature, it 
becomes 




Another important quantity is the Wigner delay time 
which gives the time delay of the scattered electron due 
to its interaction with the scattering field ( hera-the os- 
cillating potential). It relates to the S-matrix bycJ 

ME) = -f r-K^fg] = -£±H*r>. (15) 

where N c is the number of open channels. Physically, 
the Wigner time represents the time spent by a wave 
packet passing through the scattering region. The charge 



pumping is supposed to be adiabatic when r is much 
greater than the Wigner delay time t w . 

It is obvious that the net charge transfer in one cycle 
is zero for a single time-periodic barrier since = T<_. 
Then the simplest system which may induce the non- 
trivial charge pumping should include at least two bar- 
riers. As an example we consider a mesoscopic sys- 
tem with two time-periodic barriers connected through 
ideal leads to two electron reservoirs with the same 
chemical potentials /i. The potentials are described by 
Fi(a;,t) = V° + V lC os{iot), V 2 {x,t) = V§ , and V 3 (x,t) = 
V 3 ° + V3Cos(u)t + (f>). This appears to be a simplified 
model for the Switkes et aFs experiment, nevertheless 
it turns out that some essential characteristics can be 
exhibited, as we will address below. In the following nu- 
merical calculations, fi = 75 mev, m* = 0.067m e (with 
m e as the mass of the free electron), V 2 ° = — 30 mev, 
V% = 0, and N c is determined by a natural condition: 
\ T ^? + \R^\ 2 - L < c e with c e (= 1.0 x 10~ 4 in this 
paper) as a defined error. 
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FIG. 1. The pumped currents I(<f> — n/2) versus the bar- 
rier height V for different pumping frequencies. Dotted lines 
fit the relation 7(y) oc V 2 . 

The general characteristics of quantum pumping are 
shown in Figs.l and 2. The parameters in Figs.l, 2, and 
3 are chosen as V° = V 3 ° = 50 mev, and Vi = V3 = V. 
Figure 1 shows that the pumped current I(V) is propor- 
tional to V 2 for small pumping amplitude V, with the 
proportional factor depending on the driving frequency 



3 



(ft = 1). But it deviates from ^-dependence for the 
strong pumping case. On the other hand, the pumped 
current is sinusoidal dependence on <f> for weak pump- 
ing, and becomes nonsinusoidal dependence on <f> when 
V increases, as seen in Fig. 2. Another important char- 
acteristic shown in Fig. 2 is that I(ir + (f>, V) = —I(<j>, V) 
for all amplitude strengths, and \I((j>, V)\ is maximum 
at <j) = 7r/2 or <j) — 3n/2 for weak pumping. Remark- 
ably, these results for /(</>, V) are in agreement with the 
experimental observation in Ref. [3] . 



V.7.0 

V.11.0 

V=14.5 




ferent from that of the incoming wave. Consequently, 
the adiabatic condition, which requires that the quantum 
state is at the same instant state in the whole evolution, is. 
not satisfied. Note that the formula derived by Brouwem 
may be valid merely under the adiabatic condition, and 
thus the method developed here may be quite useful. 




Fig.3 
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FIG. 3. The pumped current and the Wigner delay time 
versus the insert energy Eq for V = 7.0 and u = 6.0. 



FIG. 2. The pumped currents versus the phase difference 
<j> for three different V and ui = 3.0. 

Figure 3 shows that sharp peaks in the Wigner time oc- 
cur at the resonance insert energies Eq = rihto. Besides, 
jumps in the pumped current as a function of E appear 
at the peak positions of the Wigner time. The direction 
of the current depends crucially on the insert energy. It 
is interesting to note that the adiabatic condition is not 
necessary in our calculations. Figure 3 indicates that the 
maximum value of t w is about 5.5 ns for Tilu = 6.0 mev 
(corresponding to t ~ 0.7 ps ), which is much greater 
than the pumping cyclic time r. Then we may say that 
the method described here is beyond adiabaticity. Ac- 
tually, the nonadiabatic effects are only important for 
the strongly photon-assisted transport since t w is greater 
than r only if the energy of the incoming wave is approxi- 
mately equal to the resonance energy for photon-assisted 
tunneling. Physically, by emitting or absorbing photons, 
the outgoing waves may be at the quantum states dif- 



It is quilt intriguing to note from Fig. 4 that I{4> = 
0) is nonzero for V\ ^ Va, while the corresponding 
areas enclosed by the path in the parameter space 
{Vi(a;, t), V%{x, t)} are zero. Although the pumped cur- 
rents in the above case were predicated to be zero under 
the adiabatic approximation, the deviation from zero is 
reported experimentally at strong pumping,tl just as we 
observed here in terms of a rigorous theoretical analy- 
sis which is also valid for strong pumping. Moreover, 
from comparison with that I(4> = 0) = for V\ = V3, it 
is now clear that the present nonzero pumped currents 
stem from the spatial asymmetry of potentials V\ =/= V3, 
which is-.coincident with the result obtained by Wagner 
in Ref.tll: the nonzero currents may be observed in a 
single osscillating potential but with asymmetric static 
potential. Actually, to observe a pump current at zero 
applied bias, it seems that the inversion symmetry should 
be broken, either in real or in k space. 



4 




Fig.4 



FIG. 4. The pumped currents versus the phase difference 
<j> for lj = 3.0 for different Vi and V3. 

The fact that I((f> = 0) is nonzero at strong pumping 
may be understood based on a scenario of the nonadi- 
abatic geometric phase.li3 Pumped currents are deter- 
mined-, hy geometric phase accumulated in the evolu- 
tionMm Under the adiabatic approximation, I(cf> = 0) = 
is predicted theoretically because the corresponding 
adiabatic geometric phase is zero. While it is now clear 
that the nonadiabatic geometric phase may be nonzero 
even in the case where the area enclosed by the path 
in the parameter space is zero (thus the adiabatic phase 



is zero).t3 Therefore, the nonadiabatic correction to the 
currents should be taken into account for strong pump- 
ing whenever the adiabatic condition is not well satisfied. 
Physically, it is reasonable to believe that the observed 
nonzero pumping at phase = for the strong pump- 
ing stems from the nonadiabatic correction when the in- 
version symmetry is broken. Practically, the asymmet- 
ric spatial potential might be present in the experiment, 
which may originate from either the shape-distorting ac 
voltages, eg from the internal potential established during 
transport .0 Since the current calculated in this approach 
is conserved since \T^±\ 2 + l-fi^l 2 = 1.0, no internal po- 
tential appears explicitly in the present formulism. It 
is worth pointing out that nonzero pumped currents for 
4> = tt are also seen in Fig.4, which seems to contradict 
with that in Refi. Also note that a nonzero I{4> = 7r) 
was also predicted by another totally different theoreti- 
cal studyO, so this contradication is still an interesting 



open question at present. 

In summary, we developed a method to calculate the 
pumped current and Wigner delay time in a mesoscopic 
system with a series of time-periodic barriers connected 
to two electron reservoirs, which appears to be applicable 
for strong pumping cases. 
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the University of Hong Kong. 
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